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ABSTRACT 

Type la supernovae are generally thought to be due to the thermonuclear explosions of carbon-oxygen white 
dwarfs with masses near the Chandrasekhar mass. This scenario, however, has two long-standing problems. 
First, the explosions do not naturally produce the correct mix of elements, but have to be finely tuned to pro- 
ceed from sub-sonic deflagration to super-sonic detonation. Second, population models and observations give 
formation rates of near-Chandrasekhar white dwarfs that are far too small. Here, we suggest that type la su- 
pernovae instead result from mergers of roughly equal-mass carbon-oxygen white dwarfs, including those that 
produce sub-Chandrasekhar mass remnants. Numerical studies of such mergers have shown that the remnants 
consist of rapidly rotating cores that contain most of the mass and are hottest in the center, surrounded by dense, 
small disks. We argue that the disks accrete quickly, and that the resulting compressional heating likely leads 
to central carbon ignition. This ignition occurs at densities for which pure detonations lead to events similar to 
type la supernovae. With this merger scenario, we can understand the type la rates, and have plausible reasons 
for the observed range in luminosity and for the bias of more luminous supernovae towards younger popula- 
tions. We speculate that explosions of white dwarfs slowly brought to the Chandrasekhar limit — which should 
also occur — are responsible for some of the "atypical" type la supernovae. 
Subject headings: binaries: close — supernovae: general — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SN la) result from thermo-nuclear ex- 
plosions of carbon-oxygen white dwarfs (CO WDs). They are 
generally thought to be triggered when the WD approaches 
(for accretion) or exceeds (for a merger) the Chandrasekhar 
mass, and the density and temperature become high enough 
to start runaway carbon fusion. This scenario, however, nei- 
ther naturally leads to explosions that reproduce the observed 
lightcurves and remnants, nor easily accounts for the varia- 
tion in SN la properties and their dependence on host galaxy. 
Furthermore, the predicted formation rates are lower than ob- 
served. 

The above leads us to reconsider the assumptions under- 
lying the standard picture. After reviewing the salient prop- 
erties of SN la (Sect. |2j, we first argue that their rates are 
easiest to understand if most mergers of CO WDs lead to 
SN la, independent of whether or not the total mass exceeds 
the Chandrasekhar mass. After a brief discussion of previous 
sub-Chandrasekhar models (Sect. @), we next argue that igni- 
tion following mergers is likely (Sect.|5]l. We close with some 
ramifications (Sect. [6}. 

2. PROPERTIES OF INDIVIDUAL SN IA AND THEIR REMNANTS 

The lightcurves of most SN la are remarkably similar, and 
can be described empirically as a (nearly) single-parameter 
family, in which timescale and maximum luminosity are 
tightly correlated, with longer-lasting explosions being more 
luminous and energetic (Phillips 1993). Underlying this vari- 
ation is the amount of radioactive 56 Ni. From SN la spec- 
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tra, iMazzali etaT] d2007l) find that this ranges from — 0.1 to 
0.9M Q . They also infer ^0.1M Q of stable iron-peak el- 
ements, and an amount of intermediate-mass elements that 
is such that the total mass of nu clear processed materi al is 
roughly constant, just over 1M . Stritz inger et alj (120061) use 
peak luminosities to infer similar 56 Ni masses, but their total 
masses, inferred from the times that the ejecta become opti- 
cally thin, do not cluster, but range from 0.5 to 1.3M . 

From SN la spectra, it is also clear that the ejecta are 
stratified, with iron -peak elements formed deeper inside 
(IMazzali et al.l 120071) . a nd (small amounts of) unprocessed 
carbon on the outside ((Thomas et al. 2007). Hydrogen is 
absent (Leonard 2007|), inconsistent with expectations for a 
hydrogen-rich progenitor companion with a strong wind or 
easily entrained envelope. 

Studi es of SN la remnan ts paint a similar picture. For in- 
stance, Ba denes et alj d2006l) find that for Tycho, models with 
about twice as much iron-peak as intermediate-mass elements 
best reproduce the X-ray spectrum, consistent with SN 1572 A 
having been a "standard" SN la (confirmed beautifully using 
light echoes: iRrause et al.ll2.008h . From the ionization struc- 
ture, iBadenef^ninHQQ^) infer stratified ejecta, strongly sug- 
gesting the explosion was (partly) supersonic. Comparing 
the predicted remnant flux, Chandr asekhar-mass models are 
too luminous. Bade nes~et alj |2006) attribute this to a break- 
down in their one-dimensional remnant models, but it could 
also indicate a lower mass. A separate clue is that most rem- 
nants appear to be evolving into a constant-density medium, 
with properties like those of the warm interstellar phase, and 
unlike those expecte d for progenitors with strong, fast winds 
(Bad enes etai1l2007l) . 

The rates and properties of SN la depend on environ- 
ment, with star-forming galaxies having higher rates of, on 
avera ge, brighter SN la than passively evolving galaxies 
(e.g.jHamuv et alJl995UMannucci et aLlbooBllSullivan et alj 
120101) . The rate has bee n suggested to depend on both mass 
and star-formation rate (Man nucci et alj 120051) . or to simply 
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be a roughly consta nt fraction, of — 1%, of the in stantaneous 
WD fo rmation rate dPritchet et al .1120081: but see Ma oz et alj 
l2010al whose data suggest a bre ak beyond — 2Gyr [their 
Fig. 5]). Consistent with the latter. Ra skin et al.1 (12009) found 
that SN la are delayed by -200 to 500 Myr from the onset 
of star formation (as inferred from local environments). In- 
tegrated over a Hubble time, - 0.0023 ± 0.0006 SN la seem 
to occur for every solar mass formed (Ma nnucci et al.ll2005b 
iMaoz et al]|2010al) . with even higher numbers, of > 0.0034 
inferre d from galaxy cluster iron abundances (Maoz et al. 
l20T0bl) . 

In summary, the lightcurves, spectra, and remnants of SN la 
supernovae seem to require center-lit explosions of CO WDs, 
with masses of > IMq. Their progenitor systems likely 
did not host companions with strong hydrogen-rich winds or 
loosely bound envelopes. The rate of SN la appears to be 
— 1 % of the WD formation rate, and younger systems pro- 
duce more luminous explosions. 

3. EXPECTED SN IA RATES 

The number of SN la per solar mass formed, —0.0023, is 
higher than expected for Chandrasekhar-mass systems, both 
from counts of suitable intermediate-mass progenitors (e.g., 
lMaozj|2008h and from po pulation synthesi s calculations (for 
recen t work, see, e.g., iRuiter et"al] 120091: iMenneke ns et alj 
2010). The details are complex and metallicity dependent, but 
below we elucidate the issues with rates e stimated from b asic 
principles (for a more formal analysis, see|Greggio 20051). 

For the primary in an interacting binary to leave a CO WD, 
it must be massive enough not to leave a He WD, but not 
so massive that it forms an ONe WD or neutron star. Using 
Figs. 1 and 2 of Webbinkj (120081) . we estimate a mass range 
1 .8 < Mi < 7M Q (the lower limit exceeds that for single stars 
b ecause the star n eeds to ignite helium after its interaction). 
A lChabrierl d2005l) initial mass function produces ni^<M,<7 = 
0.067 such stars per solar mass formed. To produce a CO WD 
of >0.7M Q (half the Chandrasekhar mass), a >3.5M© pri- 
mary is required, of which only «3.s<m,<7 = 0.020 are formed. 
For these suitably massive stars, a fraction fc m ~ | is in bina- 
ries, and, of those, a fr action /io<p<2ooo — 0.2 have periods be- 
tween 10 and 2000 d dDuquennov & Mavoril 19911) . such that 
they interact (P < 2000 d), but only after the main sequence, 
with a fully formed helium core (P > lOd). 

The further evolution depends on whether mass transfer is 
stable or not. If it is unstable, common-envelope evolution 
will drastically shrink the orbit, leading to possible further 
evolu tion via the single-degenerate channel ( W helan & Ibenl 
1973'). If it is stable, the system remains wide and further evo- 
lution leads to a second CO WD as well as a second, unstable 
mass-transfer phase that shri nks the orbit, as required for the 
double-degenerate scenario (Webbink 1984; Iben & Tutukovl 

cm. 

For the giants considered here, mass transfer to a less mas- 
sive compa nion is generally expected to be dynamically un- 
stable (e.g., Webbink 2008). So one naively expects a small 
fraction / w id e left in wide orbits, and a near-unity fraction 
/dose left in close orbits. Empirically, howeve r, the first mass- 
transf er phase sometimes le aves wide orbits (Nelemans et al. 
2000; for a discussion, see Webbink 2008), presumably for 
nearly equal-mass binaries (w hich may be relatively common; 
Pinsonneault & Stanek 2006). Indeed, the existence of fair 
numbers of both double-degenerate binaries and cataclysmic 
variables suggests neither fraction is small. Below, we as- 
sume / w ide — 5 and / c i ose ~ §; likely, neither is off by more 



than 50%. 

3.1. Single degenerates 

In principle, the number of CO WDs formed 
in close orbits with non-degenerate companions, 
«i.8<M,<7/bin/io<p<200o/ciose ^ 0.006 per solar mass, could 
reproduce the S N la rate, and many rou tes to explosions have 
been proposed ( Iben & Tutukov 1984). No route, however, 
seems both common and efficient in growing the WD to the 
Chandrasekhar mass. If mass transfer is too slow, novae 
occur, which appear to remo ve as much mass as was accreted 
dTownslev & Bildsten] l2004t possible counterexamples are 
RS Oph and U Sco). If it is faster, hydrogen burns stably, 
but only a small r ange avoids expansion and mass loss 
(iNomoto et al]l2007l) . 

Empirically, t he only efficient system s appear to be the su- 
pers oft sources (Rappap ort et al.lll994|). but thos e are far too 
rare dDi Stefanol2010l;lGilfanov & Bogdanl2010l) . We may be 
missing systems, e.g., rapidly accreting WDs that expanded 
and hid from X-ray view. However, for such sources — as for 
many single-degenerate channels — the absence of evidence 
for hydrogen and wind-blown bubbles is surprising. The lack 
of convincing solutions to these issues motivates us to look 
for alternative progenitors. 

3.2. Double degenerates 

Given the near-unity mass ratio required to keep a wide or- 
bit in the first mass-transfer phase, the two CO WDs are ex- 
pected to have similar masses. To estimate the number of dou- 
ble degenerates with total mass exceeding the Chandrasekhar 
mass, we thus use the number of binaries with sufficiently 
massive primaries, «3.5<M,<7./bin/io<.p<200o/wide ^ 0.0009 per 
solar mass formed. This is less than half the required number, 
which poses a significant problem, especially as some systems 
will be too wide to merge in a Hubble time. 

Indeed, this realisation motivates our consideration of sub- 
Chandrasekhar merger models: if all mergers of CO WDs 
would lead to SN la, one has ni.8<M,<7/bin/io</><2ooo/wide ^ 
0.003 possible progenitors, which is consistent with the ob- 
servations. 

Furthermore, sub-Chandrasekhar mergers could explain the 
observed delay-time distribution. Generally, distributions of 
formation times are shallower, oc ?~ 5 (Pritchet et al. 2008), 
than those of merger times, oc f" 1 (e.g., Greggio 2005). Thus, 
one expects the SN la rate to scale with the WD formation 
rate. Quantitatively, the scale factor is —0.01 dPritchet et alj 
2008), while the fraction of WDs formed in double degener- 
ates is about /bin/io</><2ooo/wide ^ 0.04. Thus, - 25% of the 
double degenerates should merge fast compared to the pro- 
genitor lifetime. But the scaling with WD formation rate will 
hold only for a duration roughly equal to the lifetime of the 
lowest-mass progenitor. For >3.5M Q progenitors, this will 
be —200 Myr, much shorter than observed, but for > 1.8M 
stars, it is — 1.7Gyr, consistent with the observations. 

4. PREVIOUS SUB-CHANDRASEKHAR MODELS 

We are n ot the first to consider su b-Chandrasekhar models 
for SN la. Woosley & W eaver! d 19941) suggested that single- 
degenerate sub-Chandrasekhar mass explosions might be trig- 
gered by detonations of overlying He layers. However, those 
lead to stratifications inconsistent with the observations (see 
Sect.©. 

In contrast, ISim et all d2010l) studied central detona- 
tions of sub-Chandrasekhar WDs (ignoring the ignition 
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Figure 1. Temperate-density profile of the remnant of a merger of two 
O.6M0 C O WDs. The thin cyan cu rve shows the profile from the sim- 
ulations of Loren-Aguilar et al. (2009), and the thick black curve the one 
that results from mixing the central convectively unstable region. After its 
formation, the remnant will accrete and become denser and hotter, evolv- 
ing at constant entropy (green dashed, diagonal contours) as long as the 
accretion timescale is shorter than t u =CyT /e u and tqq = CyT /eqq, the 
timescales on which neutrinos can cool and carbon fusion can heat the core 
(blue dotted and red dashed contours, respectively). The magenta solid curve 
shows where ignition occurs (t„ = tqq) and the thick black dot-dashed one 
where degeneracy is lifted (normalized electron chemical potential r\ = 0). 
(The entropies S, specific heats CV, and gain/loss rates e were calculated 
using routines from the stellar evolution code MESA |Paxtonetal. 2010; 
http://mesa.sourceforge.net |; uncertainties in ecc. while large, do not influ- 
ence our conclusions) 

mechanism, which we address in Sect. [5j. They found 
lightcurves and spectra similar to SN la, and showed that 
for more massive WDs, more iron-peak elements were cre- 
ated, with M Ni ~ {0.06,0.32,0.58, 0.85}M o for M WD = 
{0.88,0.97, 1.06, 1.15}M Q . 

That mass dependence arises because iron-peak elements 
are produced only in regions that reach ~4 x 10 9 K be- 
fore degeneracy is lifted, which requires a density p > 
10 7 gcm~ 3 (to produce intermediate mass elements requires 
p>2x 10 6 gcm" 3 ). For the same reason, detonations of near- 
Chandrasekhar WDs produce too much nickel: almost the 
whole WD is above the critical density. This conundrum can 
only be solved by a n initial deflagra tion phase, which allows 
the WD to expand (Kho kMovlll991l) . For lower mass WDs, 
this is not necessary. 

5. IGNITION 

We argued that the SN la rate and delay-time distribution 
could be understood if mergers of CO WDs lead to SN la even 
for sub-Chandrasekhar total mass. If mergers lead to explo- 
sions similar to the detonat ions of sub-Cha ndrasekhar WDs, 
they will appear like SN la dSim et alj|2010l) . Assuming more 
massive mergers also produce more 56 Ni, the range in lumi- 
nosity and duration follows naturally, as does the correlation 
with parent population age. 

The critical remaining question is whether merger products 
become sufficiently hot to ignite. This has been addressed 
partly by recent simulations, which include caref ul treatment 
of the equation of state and of nuclear burning (|Yoon et ail 
120071: lLoren-Aguilar et aT]|2009t IPakmor et alJlfoiOl) These 



show qualitative differences between mergers of unequal and 
equal mass binaries^ If one WD is significantly lighter (and 
thus larger), mass transfer leads to its total disruption, with 
the material wrapped around the more massive companion; 
the merger product has a core that is cooler and rotates more 
slowly than the envelope. Any nuclear processing happened 
at the co re-envelope interface, a s likely would any subsequent 
ignition dNomoto & Ibenlll985b lYoon et al J 120071) . Thus, we 
do not think such mergers lead to SN la. 

If instead the two WDs have more equal masses, as ex- 
pected for CO+CO binaries (Sect. 13. 2\ , the merger remnant 
is fully mixed and hottest in the center. Initially, it rotates 
differentially, but this is dissipated, and one is left with a 
core holding ~ 80% of the mass and rotating at a uniform 
rate near the mass-shedding limit, surrounded by a somewhat 
sub-Keplerian, very dense, partially degeneracy-pressure sup- 
ported "disk" with a steep surface density gradie n t (E o c r~ 5 ). 

In the simulation s of lLoren-Aguilar et ail (12009b and 
Pakm or et al.l (120101) . the mergers of equal-mass WDs do 
not become hot enough to ignite carbon burning, except for 
masses above ~ 0.9M Q . The explosion of the resulting 
> 1.8Mq remnant, however, leads to a subluminous SN la. 
This is not surprising: generally, merger remnants have cen- 
tral densities similar to those of the (more massive of the) 
pre-merger WDs. For a 0.9M Q WD, the central density is 
~ 1.5 x 10 7 gcm~ 3 , and thus little 56 Ni will be produced (see 
Sect. 0J, leading to a subluminous explosion^ 

We now turn to the simulation o f the m erger of two 
0.6M Q WDs of lLoren-Aguilar et all d2009l) . This mass 
is close to the empirical m ean mass of WDs (^0.65M Q ; 
iTremblav & Bergeronl 120091) . and thus, in our picture, the 
merger should lead to a typical SN la. From the simula- 
tion, the remnant's central density and temperature are ~ 
2.5 x 10 6 gcm" 3 and ~6 x 10 8 K, respectively. This would 
appear to be close to what is required for ignition, but since 
the central ^O.O7M is convectively unstable (Figure[TJ, the 
temperature will be reduced rapidly, on a convective turn- 
over time. The region's average entropy, s/k ~ 16 per ion 
and ~ 1.3 per electron, corresponds to a core temperature of 
-3.8xl0 8 K. 

While this merger remnant is too cold (and insufficiently 
dense) to produce a SN la, it could explode later. One pos- 
sible heating mechanism is the accretion of the thick disk, 
which will lead to an increase in density and thus to com- 
pressional heating. Assuming the source remains at the mass- 
shedding limit and in roug hly solid-body rotation (e.g., due to 
magnetic fields; Sect. l5.1l below). the density will increase to 
~3 x 10 7 gcm" 3 (for a cold WD: iGerovannis & H adiopoulos 
1989). If the contraction is fast enough to be (nearly) adia- 
batic, the central temperature would increase to ~1.4x 10 9 K 
(Figure[TJ, easily hot enough to ignite. 

The timescale for compressional heating is th = 
(M core /M disk )r acc , where M core /M disk ~ 5 is the core-to- 
disk mass ratio, and T acc is the disk accretion timescale, which 
we estimate using the usual a formalism, 

_ A/disk _ _!/r disk \2 



5 What constitutes "equal" is not yet known, though 0.8 and 0.6Mq is 
not. For reference, in the range 0.6-0.9Mq, to have central densities within 
afactor of two requires AM <Q.13Mq. 

6 The rare mergers of even more massive WDs should lead to more lumi- 
nous explosions. 
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Here, M is the accretion rate, and rdi s k — O.O27? , h ~ 
O.OO55/?0, and Td yn = ~ 5 s are the disk radius, scale 
height, and dynamical ti me (with numerical values from 
lLoren-Aguilar et al] 120091; note that their Table 1 lists outer 
disk radii; their Figure 3 shows that most mass is at much 
smaller radius). 

Thus, we find that the remnant is heated on a timescale 
th ~ 5r acc ~ lOh, which is much shorter than the neutrino 
cooling and carbon burning timescales even at the ignition 
line (where both are ^3 x 10 3 yr, see Figure [TJ. Therefore, 
neutrino cooling can be ignored, and compressional heat- 
ing will continue until the disk is exhausted or the fusion 
timescale has become shorter than the accretion timescale. 
From Figure Q] the latter happens when p ~ 1.6 x 10 7 gcm~ 3 
and T ~ 10 9 K, which is slightly before disk exhaustion (at 
p ~ 3 x 10 7 gcm" 3 ; see above). At this point, a nuclear run- 
away is inevitable. 

5.1. Complications 

Above, we argued it is plausible that merger remnants will 
heat up sufficiently to ignite carbon, but we made a number 
of simplifying assumptions that deserve further study. First, 
the "alpha" formalism may be inappropriate for estimating the 
accretion timescale for a small, massive, and thick disk. In- 
stead, the relevant timescale m ay be the much lo nger cool- 
ing timescale of the envelope (lYoon et al.l 12007!) . If trans- 
port of angular momentum is the determining factor, how- 
ever, our timescale is the correct order-of-magnitude estimate. 
Second, as the accretion rate is highly super-Eddington, a 
(strong) wind may form, which may diminish or enhance the 
compression depending on its specific angular momentum. 
Third, accretion could heat the envelope significantly. For 
cold WDs, rapid acc retion leads to off-center ignition (e.g., 
iNomoto & iben|[T985l) . but it likely is less important when the 
WD core is hot. Fourth, we assumed the remnant core ro- 
tates r oughly uniformly. This is found in the merger simula- 
tions dLoren-Ag uilar et aTll2009[ and references therein), but 
may reflect artificial viscosity associated with smooth parti- 
cle hydrodynamics. A strongly differentially rotating remnant 
would be less dense and suffer less from compressional heat- 
ing. 

If differential rotation is present, it also leads to additional 
effects that we ignored. It drives a number of processes 
that tend to eliminate it (Piro 2008). In particular, it will 
wind up magnetic fields until their energy is of order the 
differential rotation energy, B w (mAft/tf 3 ) 1 / 2 . With IQ ~ 
10 50 gcm 2 s _I and R ~ 0.0125 Rq, the inferred field ranges 
from ~ 10 9 G if differential rotation is driven by the accretion 
(AQ ~ TjJ 1 ~ 10" 5 s _1 ) to ~ 10 n G if it is due to the merger 
(A57 ~ ft ~ 10 _1 s" 1 ). Empirical evidence for field generation 
comes from arguably the best candidate WD merger remnant, 
RE J03 17-853. This WD is massive, M ~ 1.35M Q E spins 
rapidl y, P = 725 s, and h as a strong, B ~ 340 MG magnetic 
field (iBarstow et al.l fl995). Such strong fields, if they emerge 
sufficiently fast, would also spin-down the WD through mag- 
netic dipole emission or through coupling with the accretion 
disk and/or a wind (as may have happened for RE 03 17-853). 

7 In the context of our scenario, the existence of this object is puzzling. It 
may be the result of an unequal mass merger. 



The concomitant loss of rotational support would add to, or 
might even dominate, the compressional heating due to accre- 
tion. 

A separate issue is that to produce SN la, we require det- 
onatio ns. Spontaneous detonations are very difficult i nitiate 
(e.g JNiemeyer & Woosleylll997l:ISeitenzahl et alJl2009h . al- 
though our case is helped by having a large super-critical re- 
gion. Furthermore, at our relatively low central density, burn- 
ing is in the distributed rather than flamelet regime, which 
may help a possible d eflagration transition to a detonation 
(Nie meyer &Woosleyll 19971) . 

A final complication is the assumed composition. Merger 
simulations so far have used carbon and oxygen only, but 
~ 1% of the mass will be helium. If this is burned during 
the merger, about 10 49 erg would be generated, roughly dou- 
bling the thermal energy of the remnant and possibly leading 
to carbon ignition. 

6. CONCLUSIONS 

We have argued that type la supernovae result gener- 
ally from mergers of CO WDs, even those with sub- 
Chandrasekhar total mass. If true, a number of interesting 
consequences arise. First, the merging WDs should have to- 
tal masses between double the lowest and highest possible 
CO WD masses, i.e., 1 < M tot < 2.4M (though ignition dur- 
ing the merger may cause a break in properties at > 1.8M Q ; 
Pakm oretal.1 120101) . This range could account for "super- 
Chan drasekhar" SN la (e.g.. lHowell et alj|2006t IScalzo et alj 
2010). As the typical total mass will depend the population's 
age, it also explains the empirical age-luminosity relation. 

Second, the rate of mergers — and thus of SN la — may 
show a break at ~ 1.7Gyr (the lifetime of a ~ 1.8M star), 
likely being more sensitiv e to the WD formation rate before- 
hand dPritchet et alj20 08). and to the merger-time distribution 
thereaf ter. Suc h a break may have been observed (Ma oz et ail 
2010a; but see lMaoz et aTll2010bl) . 

Third, ignition will likely take place in rapidly rotating 
objects, (possibly) surrounded by small disks. This may 
have interesting con sequences for the explosion dynamics 
dPfannes et al J 12010). the initial shock breako ut dPiro et alj 
120101) . early-time spectra dMazzali etafll2005T) . and the su- 
pernova remnants. 

Finally, we speculate that WDs slowly pushed to the 
Chandrasekhar mass — such as should be produced in single- 
degenerate systems — are partly responsible for the population 
of "atypical" SN la. Further contributions to that popula- 
tion might come from mergers that ignite during the merger 
proper, and from unequal mass mergers that ignite off-center 
and/or explode only partially. 
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Berro and Pablo Loren-Aguilar for providing details on their 
merger simulations. We are grateful to Bill Paxton for making 
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